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A search for the production of neutral Higgs bosons decaying into r r~ final states is pre- 
sented. One of the two r leptons is required to decay into a muon. The data were collected 
by the D0 detector and correspond to an integrated luminosity of about 1.0 fb~^. No excess 
is observed above the expected backgrounds. The results are interpreted in the Minimal Su- 
persymmetric Standard Model. In the mass range 90 < raA < 200 GeV values of tan /3 larger 
than 40-60 are excluded for the no-mixing and the m™^" benchmark scenarios. 

1 Introduction 

The contribution of r^r^ final states from the Standard Model (SM) Higgs production is too 
small to play any role in the SM Higgs searches in pp collision at the Tevatron due to the 
large irreducible background from Z t^t" production. This is different in the Minimal 
Supersymmetric Standard Model (MSSM), which introduces two Higgs doublets leading to five 
Higgs bosons: a pair of charged Higgs boson (H^); two neutral CP-even Higgs bosons (h,H) and 
a CP-odd Higgs boson (A). At tree level, the Higgs sector of the MSSM is fully described by 
two parameters, which are chosen to be the mass of the CP-odd Higgs, rriA, and the ratio of the 
vacuum expectation values of the two Higgs doublets, tan/3. The Higgs production cross-section 
is enhanced in the region of high tan . In the low niA, high tan (3 region of the parameter space, 
Tevatron searches can therefore probe several MSSM benchmark scenarios extending the search 



regions covered by LEP ^. Inclusive searches for (p{= H, h, A) tt have been performed with 
integrated luminosities of L = 350 pb^^ by D0 and L = 310 pb^^ by CDF^. Both searches 
require at least one r lepton to decay into an electron (rg) or a muon {ti^). In this analysis, 
only the decay r^r is considered using an integrated luminosity of L = 1.0 fb^^. CDF 
has also recently released a preliminary result with L = 1.0 fb~^^. The biggest improvement in 
sensitivity compared to the previous analysis comes from using a neural network to improve the 
separation between signal and background. 

2 Event Preselection 

The preselection requires one isolated muon with, pf^ > 15 GeV. The event is required to have 
no other muon that is matched to a track in the central detector with pi^ > 10 GeV. 

Hadronically decaying taus are characterized by a narrow isolated jet that is associated with 
three or less tracks. Three types of hadronically decaying taus are distinguished: 

Type 1: Calorimeter energy cluster, with one associated track and no electromagnetic sub-cluster. 
This corresponds mainly to the decay tt^u. 

Type 2: Calorimeter energy cluster, with one associated track and at least one electromagnetic 
sub-cluster. This corresponds mainly to the decay — > tt^tt^v. 

Type 3: Calorimeter energy cluster, with three associated tracks, with an invariant mass below 
1.7 GeV. This corresponds mainly to the decays 'k^'k^'it^{'k^)v. 

Tau decays into electrons are usually reconstructed as type-2 taus. These are not removed from 
the sample. The event is required to contain a r candidate at a distance Ai? > 0.5 from the muon 

direction. The charge of the r candidate must be opposite to the muon charge. The transverse 
momentum pj. of the r candidate measured in the calorimeter must be greater than 15 GeV 
for r-type 1 and 2, and greater than 20 GeV for r-type 3. At the same time the transverse 
momentum of the track associated with the r candidate is required to be pT > 15 GeV for 
r-type 1 and px > 5 GeV for r-type 2. In the case of r-type 3, the scalar sum of the transverse 
momenta of all associated tracks must be greater than 15 GeV. 

3 W — > fj,i/ and Multi-jet Background Estimation 

The shape of the W /iz^ + jet background distribution, where the jet is misidentified as a tau, 
was simulated using PYTHIA. The normalization, however, was obtained using data. 

A contribution to the background is expected from heavy flavour multi-jet events, where a 
muon from a semi-leptonic decay passes the isolation requirement and a jet is mis-identified as a 
r candidate. In addition, a contribution is expected from light quark multi-jet events where the 
jets fake both the tau and the muon. The multi-jet background shape is taken from events with at 
least one muon and one r candidate where the muon failed the calorimeter isolation requirement. 
The normalization of this semi-isolated sample was obtained in a multi-jet enriched sample. 

4 Final Event Selection 

A set of neural networks, one for each tau type, has been developed to separate the tau leptons 
from jets. These neural networks make use of input variables that exploit the tau signature such 
as longitudinal and transverse shower shapes and isolation in the calorimeter and the tracker. 
The neural network is trained using tau MC events as signal and multi-jet events from data as 
background to produce a variable that peaks near one for real taus and zero for jets. The tau 



candidate is required to have a neural network output greater than 0.9. In the case of type-3 
taus this is tightened to 0.95 due to the larger multi-jet background. 

It is also possible for muons to fake type-1 or type-2 tau candidates. These fakes are removed 
by ensuring that type-1 or type-2 tau candidates do not match to a reconstructed muon within 
a cone of radius Ai?^,- = 0.5. 

After selecting events with a high neural network output, there is still a considerable amount 
of background from W -|- jets production. To remove these events, the reconstructed W boson 
mass, Mw = \/2E'^ Et^{l — cos A0) is used, where = }^p>^/pj. is the estimated neutrino 
energy, calculated using the ratio of the muon momentum and muon transverse momentum 
pt^. For real W boson decays, this variable peaks near the W boson mass, whereas for the signal 
and the Z — > rr background the variable peaks at zero. Events with Mw > 20 GeV are rejected. 

To achieve the best separation of the signal from background, neural networks were trained 
for different Higgs mass points using kinematical variables. The distribution of the visible mass 
Mvis and the optimised neural networks for a Higgs mass of 160 GeV is shown in Figure 1. There 
is good agreement between the background expectation and the data. 




Figure 1: Distribution of a) the visible mass Mvis and b) neural network output distribution for a Higgs mass of 
160 GeV with all selections applied. The data, shown with error bars, are compared to the sum of the expected 
backgrounds. Overflow events are added to the last bin. Also shown, in light green, is the signal for a Higgs mass 
of 160 GeV, normalized to a cross-section of 10 pb. The systematic uncertainty on the background normalisation 

is 10% and is shown by the shaded area. 

5 Results and Conclusion 

Limits on the cross-section for Higgs boson production times the branching fraction into tau 
leptons are derived at 95% Confidence Level (CL). The output from the neural networks, shown 
in Figure 1 for one mass point, are used in the limit calculation. The distributions for the three 
tau types are used separately. The cross-section limits are calculated with the CLg method^. 

There are various sources of systematic uncertainties that affect signal and background. The 
most important are the uncertainty on the integrated luminosity (6.1%), the trigger efficiency 
(3%), the tau energy scale (1 — 11%), the uncertainty in the signal acceptance due to choice of 
parton distribution function (3.9 — 4.6%), the uncertainty of the tau track matching efficiency 
(4%), the uncertainty on the tau reconstruction efficiency (3%), the theoretical uncertainty on 
the Z cross-section (5%) and the uncertainty on the modeling of the multi-jet background (3%). 
All systematic uncertainties are included in the calculation of the expected and observed limits, 
assuming 100% correlation between signal and background where appropriate. The expected 
and observed limits are shown in Figure 2 as a function of the hypothetical Higgs mass. 

In the MSSM, the masses and couplings of the Higgs bosons depend on tan/3 and tjia at 
tree level. Radiative corrections introduce additional dependencies on SUSY parameters. In this 



analysis, the m™'^^ and no-mixing scenarios are studied. The corresponding excluded regions in 
the tan/3 — niA plane are shown in Figure 3. The cross-section at each tan/3 — iriA point was 
calculated using FeynHiggs 2.5.1 by adding the gg ^ <f) and hh (f) cross-sections for a given 
ruA- 

In the mass region 90 < rriA < 200 GeV, tan (3 values larger than 40-60 are excluded for the 
no-mixing and the m^^^ benchmark scenarios. These results are the most constraining limits 
from the Higgs to t~^t~ decay channel to date. 




Figure 2: Observed and expected 95% Confidence Level upper limit on the cross-section times branching ratio, 
using the neural network shown on both a log scale and a linear scale. The band represents the ilcr uncertainty 
on the expected limit. Also shown is the expected limit from the recent CDF result. 




Figure 3: Excluded region in the tan /? — m,4 plane for < in a) the m™^" scenario and b) the no-mixing 
scenario and excluded region in the tan /3 — niA plane for ^ > in c) the m™^'' scenario and d) the no-mixing 

scenario. 
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